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ORIGINAL ARTICLE 

A Nonlinear Mixed Effects Pharmacokinetic Model for 
Dapagliflozin and Dapagliflozin 3-O-glucuronide in Renal 
or Hepatic Impairment 

J-S van derWalt\Y Hong 2 , L Zhang 2 , M Pfister 2 , DW Boulton 2 and MO Karlsson 1 

Dapagliflozin is a sodium-glucose co-transporter 2 inhibitor in development for the treatment of type 2 diabetes mellitus. 
A semi-mechanistic population pharmacokinetic (PK) model was developed for dapagliflozin and its inactive metabolite 
dapagliflozin 3-O-glucuronide (D30G) with emphasis on renal and hepatic contribution to dapagliflozin metabolism. Renal and 
hepatic impairment decreased the clearance of dapagliflozin to D30G and the clearance of D30G.The fraction of D30G formed 
via the renal route decreased from 40-55% in subjects with normal renal function (creatinine clearance (CLcr) > 80 ml/min) to 
10% in subjects with severe renal insufficiency (CLcr = 13 ml/min). The model-based simulations suggested that the increase of 
systemic exposure (AUCJ of dapagliflozin and D30G was less than twofold in subjects with mild or moderate renal impairment. 
This population modeling analysis presents a useful approach to evaluate the impact of renal and hepatic function on the PK 
of dapagliflozin. 
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Dapagliflozin is a therapy in development for the treatment of 
type 2 diabetes mellitus (T2DM). Under normal conditions, 
~180g of glucose is filtered through the renal glomerulus per 
day, and virtually all of this is reabsorbed into blood by the 
proximal tubules with the major facilitator of reabsorption being 
the sodium-glucose cotransporter. 1 Dapagliflozin acts revers- 
ibly and specifically on the sodium-glucose co-transporter 2 in 
the proximal tubules of the kidney to inhibit the reabsorption of 
glucose and to promote urinary glucose excretion. 2 

Dapagliflozin is rapidly absorbed following oral administra- 
tion with an absolute oral bioavailability of 78%. 3 The systemic 
exposures of dapagliflozin increase in a dose-dependent 
manner for doses ranging from 0.1 to 500 mg. The half-life of 
dapagliflozin (12.5h) and the corresponding sustained inhi- 
bition of urinary glucose reabsorption over 24-h after dose 
at the proposed clinical dose (10mg) is amenable to once- 
daily dosing. 2 Dapagliflozin is metabolized by uridine diphos- 
phate glucuronosyltransferase (UGT) 1A9 to dapagliflozin 
3-O-glucuronide (D30G identified as M 15 in chromatogra- 
phy), which is 2600-fold less potent than the parent drug with 
regard to sodium-glucose co-transporter 2 inhibition. Follow- 
ing a single 50-mg dose of [14C]-dapagliflozin, D30G and 
dapagliflozin represented >95% of the total radioactivity in 
all plasma samples. In healthy subjects, -61% of the admin- 
istered dose of dapagliflozin is recovered in urine as D30G. 4 
Less than 2% of the administered dose of dapagliflozin is 
excreted in urine as unchanged dapagliflozin. 2 

Approximately 40% of adults with T2DM have comorbid 
chronic renal disease, and T2DM is widely regarded as the 
leading cause of end-stage renal disease. 5 Individuals with 
T2DM have a higher incidence of liver function test abnormali- 
ties than those without T2DM, and therefore, hepatic impair- 
ment (HI) is also an important consideration in this population. 6 



A pharmacokinetic (PK)/pharmacodynamic (PD) study of 
dapagliflozin in renal impairment (NCT00554450) and a PK 
study in HI have been described. 7-9 In the renal impairment 
study, systemic exposure of both dapagliflozin and D30G 
increased incrementally with decreased renal function. The 
PK of dapagliflozin in patients with T2DM and normal renal 
function are comparable to that observed in healthy subjects. 2 
Compared with healthy subjects and T2DM subjects with nor- 
mal renal function (creatinine clearance (CL cr ) > 80 ml/min), 
the dapagliflozin exposure (area under the concentration time 
curve from time zero to infinity (AUC 0 _J) increased by 28% in 
patients with T2DM and mild renal impairment (CL cr > 50 and 
CL cr < 80 ml/min), by 52% in patients with T2DM and moder- 
ate renal impairment (CL cr > 30 and CL cr < 50 ml/min), and by 
75% in patients with T2DM and severe renal impairment not 
on hemodialysis (CL cr < 30 ml/min). The categorization of vari- 
ous stages of renal impairment was trial defined; however, the 
Kidney Disease Outcomes Quality Initiative clinical practice 
guidelines for chronic kidney disease define mild, moderate, 
and severe renal impairment as glomerular filtration rate of 
60-89, 30-59, and 15-29 ml/min, respectively. 10 

Glucuronide conjugation is usually preserved in patients 
with hepatic insufficiency until the disease is severe (Child- 
Pugh Class C). In these cases, increased systemic exposure 
of dapagliflozin and decreased D30G would be expected. 
In the HI study, patients with mild, moderate, and severe 
HI showed increases in the D30G AUC 0 _„ by 6, 100, and 
29%, respectively, compared with age-, weight-, gender-, 
and smoking status-matched healthy subjects. These values 
were highly dependent on the calculated CL cr of each group. 

Furthermore, as UGT1A9 is expressed in both the kidney 
and the liver, 11 we postulated that the UGT1A9-mediated 
metabolic clearance of dapagliflozin occurs in both of these 
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organs, and that both hepatic and renal impairment may 
impact the metabolic clearance of dapagliflozin. Therefore, 
we developed a semi-mechanistic nonlinear mixed effects 
model to gain further insights into the effects of renal and HI 
on the population PK of dapagliflozin and D30G. 

RESULTS 
Base model 

The structural model that best described the plasma and 
urine data of dapagliflozin and D30G simultaneously is 
shown in Figure 1. The parameter estimates of the base 
model are provided in Table 1 . 

The medians (interquartile range) of fraction of dapa- 
gliflozin excreted in urine (FE) and fraction of dapagliflozin 
metabolized to D30G (FM M15 ) were 1.6% (1.2-2.1%) and 
60% (53-72%), respectively. In the base model, FE and 
FM M15 were reduced in subjects with severe renal impairment 
(Supplementary Figure S1 online). The CLP M15 and the vol- 
ume of distribution of D30G (V2M) were found to be corre- 
lated (r= 0.6). Overall, the single-dose base model described 
the observed data well but underpredicted the population 
D30G concentrations in plasma in subjects with CL cr < 50 ml/ 
min (data not shown). 

The estimates of the proportional components of the 
residual variability in plasma dapagliflozin and D30G were 
20.8 and 19.5%, respectively. Replicate residual variability in 
plasma accounted for 47.9% of the variability. As expected, 
residual variability in urine was larger: 54.2 and 51.9% for 
dapagliflozin and D30G, respectively. Replicate residual vari- 
ability in urine accounted for 80.4% of the variability. 

Covariate model 

Significant covariate relationships identified during the gen- 
eral additive modeling analysis and included in the stepwise 



covariate model (SCM) building are summarized in Table 2 
and Figure 1 . The parameter estimates of the single-dose 
covariate model are shown in Table 1. Renal function influ- 
enced the CLP M15 . The fraction of D30G formed via the renal 
route decreased from 46% in a typical subject with normal 
renal function (CL cr ~100ml/min) to 10% in a typical subject 
with CL cr = 13ml/min, the lowest value in the study popula- 
tion. This relationship is illustrated in Figure 2. 

The covariate model suggested that both formation 
(CLP M15 ) and elimination (CLM) of D30G were altered with 
HI. With severe HI (Child-Pugh Class C), CLP M15 decreased 
by 41% and V2M increased by 134%. Moderate or severe HI 
(Child-Pugh Class B or C) decreased CLM and the periph- 
eral volume of distribution of dapagliflozin (V3P) by 29 and 
60%, respectively. 

During the backward elimination, the effects of age on 
the renal clearance of unchanged dapagliflozin to urine 
(CLP renal ), CLM, and V3P, and the effects of gender on CLM 
and the central volume of distribution of dapagliflozin (V2P) 
were removed. The effect of age on the metabolic clear- 
ance of dapagliflozin to unmeasured metabolites (CLP other ) 
remained significant, and the CLP olher decreased by 2% for 
every year of increasing age above 54 years. 

The covariate model described the single- and multiple- 
dose data best and improved the population-predicted 
plasma D30G concentrations compared with the base model 
(data not shown). 

Final model 

The final model incorporated gender effects on CLM and 
total dapagliflozin clearance (CLP) and resulted in a statisti- 
cally significant improvements in model fit (AOFV = -15.828 
points). In females, CLP and CLM were 16.7 and 19.6% 
lower, respectively. The parameter estimates of the final 
model are provided in Table 1 . 
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Figure 1 The structural and covariate model to describe the renal and nonrenal elimination of dapagliflozin and D30G in healthy subjects, 
T2DM subjects with normal or impaired renal function, and patients with hepatic impairment. BCL Cr IBW, baseline creatinine clearance 
calculated using ideal body weight; CLM, renal clearance of D30G; BIO, bioavailability; CLP M15 , metabolic clearance of dapagliflozin to D30G; 
CLP Mher , metabolic clearance of dapagliflozin to unmeasured metabolites; CLP , renal clearance of unchanged dapagliflozin to urine; MTT, 
mean transit time; N, number of transit compartments; QP, intercompartmental clearance of dapagliflozin; T2DM, type 2 diabetes mellitus; 
V2P, central volume of distribution of dapagliflozin; V2M, central volume of distribution of D30G; V3P, peripheral volume of distribution of 
dapagliflozin. Dashed lines, a priori scaling; shaded areas, covariates selected during stepwise covariate model building; unshaded areas, 
added based on previous modeling experience. Covariates connected to compartments affect the relevant volume, those connected to 
pathways affect the relevant clearance. 
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Table 1 Parameter estimates 



Parameter estimates 




Single-dose base model 


Single-dose covariate model 


Final model 




Name (units) 


Typical value" 


RSE% b 


Typical value* 


RSE%" 


Typical value* 


RSE%" 


Structural fixed effects 














CLP M15 (l/h) 


6.67 


9.7 


6.88 


12.5 


7.54 


12.0 


CLP thor (l/h) 


3.83 


12.9 


4.85 


11.1 


5.35 


12.1 


CLP renal ((l/h)/(ml/minute))" 


0.00288 


10.9 


0.00286 


8.0 


0.00310 


9.9 


V2P (I) 


37.7 


4.6 


38.2 


4.3 


39.0 


4.4 


QP (l/h) 


6.59 


8.0 


6.91 


7.3 


7.07 


7.1 


V3P (I) 


61 


9.7 


69.8 


7.1 


71.5 


7.0 


MTT (h) 


0.479 


7.0 


0.47 


7.2 


0.475 


7.2 


NN 


5.69 


18.1 


5.47 


16.5 


5.45 


16.5 


BIO 


0.823 


NR a 


0.84 


NR" 


0.858 


NR d 


CLM ((l/h)/(ml/minute))° 


0.0716 


6.6 


0.0723 


6.1 


0.0799 


7.6 


V2M (I) 


2.55 


6.5 


2.22 


5.8 


2.26 


5.6 


Covariate fixed effects' 














CLP M15 -Baseline CL C[ IBW 


NA 


NA 


0.00475 


49.3 


0.00502 


42.6 


CLP M15 -Child-Pugh Class C 


NA 


NA 


-0.407 


42.8 


-0.422 


38.6 


CLP oth0I -age 


NA 


NA 


-0.0215 


24.0 


-0.0204 


22.1 


V2M- Child-Pugh Class C 


NA 


NA 


1.34 


16.3 


1.33 


17.3 


V3P- Child-Pugh Class B,C 


NA 


NA 


-0.596 


12.4 


-0.600 


12.6 


CLM- Child-Pugh Class B,C 


NA 


NA 


-0.292 


25.4 


-0.293 


23.3 


CLP-gender 


NA 


NA 


NA 


NA 


-0.167 


64.1 


CLM-gender 


NA 


NA 


NA 


NA 


-0.196 


38.9 


Random effects" 














"VCLP M15 


51.6 


15.0 


38.0 


17.4 


36.8 


16.5 


IIV.CLP othar 


52.6 


50.0 


32.9 


20.4 


31.3 


22.3 


"V.CLP^ 


62.1 


38.7 


63.1 


18.8 


61.9 


18.4 


IIV.V2P 


21.9 


32.3 


22.6 


20.0 


23.6 


20.5 


IIV.QP 


42.5 


39.5 


38.8 


41.1 


39.1 


35.9 


IIV.V3P 


57.4 


29.5 


46.5 


15.8 


46.3 


17.4 


IIV.MMT 


59.1 


28.0 


60.1 


20.5 


60.3 


21.0 


IIV.NN 


195 


13.2 


189 


10.9 


189 


10.6 


IIV. BIOS 


13.3 


NA" 


12.7 


NA d 


11.1 


15.2 


IIV.CLM 


33.3 


36.7 


27.9 


20.4 


24.9 


23.6 


IIV.V2M 


66.2 


12.8 


46.5 


13.3 


44.7 


15.5 


Correlation CLP M15 :V2M h 


r=0.60 


35.7 


r=0.56 


47.4 


r= 0.562 


49.9 


Residual error 1 














Prop, dapa_plasma 


0.208 


7.8 


0.207 


6.4 


0.207 


6.4 


Add, dapa_plasma (ng/ml) 


0.424 


87.5 


0.463 


31.7 


0.465 


29.9 


Prop, M15_ plasma 


0.195 


74.4 


0.195 


6.5 


0.195 


6.4 


Add, M15_plasma (ng/ml) 


0.633 


98.1 


0.579 


53.2 


0.585 


51.3 


Prop, dapa_urine 


0.542 


8.4 


0.565 


10.1 


0.564 


9.9 


Add, dapa urine (ng/ml) 


0.667 


506.7 


2.21 


290.5 


1.784 


640.4 


Prop, M15_urine 


0.519 


8.1 


0.564 


8.7 


0.567 


8.9 


Add, M15_urine (ng/ml) 1 ' 


0.02 fix 


0(0) 


0.02 fix 


0(0) 


0.02 fix 


0(0) 


Plasma, replicate 


0.479 


12.8 


0.473 


12.4 


0.475 


12.5 


Urine, replicate 


0.804 


6.1 


0.820 


6.0 


0.821 


6.2 


BIO, oral bioavailability of dapagliflozin; CLM, renal clearance of D30G; CLP, metabolic clearance of dapagliflozin to D30G; CLP otner , metabolic clearance of 
dapagliflozin to unmeasured metabolites; CLP ranal , renal clearance of unchanged dapagliflozin to urine; FOCE, first-order condition estimation method; IIV, 



interindividual variability; MTT, mean transit time; W, number of transit compartments; QP, inter-compartmental clearance of dapagliflozin; V2M, central volume 
of distribution of D30G; V2P, central volume of distribution of dapagliflozin; V3P, peripheral volume of distribution of dapagliflozin. 

"Random Effects parameter estimates are shown as %CV = Je°" — 1 . b RSE% is the relative standard error (SE as a percentage of estimate) obtained from 

Monte Carlo Importance Sampling EM Assisted by Mode a Posteriori estimation (IMPMAP) with Monte Carlo samples per individual = 10,000. °Renal clear- 
ances reported as l/h per ml/min creatinine clearance. "Not reported as this random effect was not MLLparameterized, a requirement for NONMEM7 IMPMAP 
method. "Fixed value. 'Covariates effects are identified by structural fixed effect then by covariate co fj ; sRaw co, (logit scale) with coefficient of variance calculated 
using SD = 8 WBI0 ■ (1 - 8-n, B | 0 ) ■ ^ja^ l0 . Correlation r = • cof y . 'Residual error parameters are identified by type (proportional (prop) or additive (add)) then 

type of observation (dapa plasma, dapa urine, D30G plasma, and D30G urine). 
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Table 2 Summary of baseline covariates and PK-covariate relationships included a priori, retained (/) or removed (x) during backwards deletion of the step- 
wise covariate model building (P= 0.01) 



Reference 
value of the 
covariate 



Baseline covariates 


Median (range) or category W/total 


relationship*'' 


CLP r„n,l 


CLP 

M15 


CLP other 


CLM 


V2P 


V3P 


V2M 


Study 


MB1 02007 


MB1 02027 


MB1 02029 


















Age (years) 


63 
(25-76) 


43 
(31-64) 


67 
(32-92) 


53.5 




X 


/ 


X 




X 




CL crlBW (ml/min)" 


60 
(13-116) 


103 
(58-143) 


40 
(20-80) 


80.14 


a priori 


/° 




a priori 








Body weight (kg) 


81.2 
(51.8-140.4) 


86.3 
(58.0- 
126.0) 


92.0 
(55.4-148.3) 


70« 


X 


a priori 


a priori 




a priori 


a priori 


a priori 


Body mass index 


19/40 


11/24 


108/163 


N/A 
















(kg/m 2 ) >30 
























Male 


24/40 


16/24 


107/163 


Male 








X 


X 






Hepatic impairment 0 
























None/mild 


40/40 


12/24" 


N/A 


None/mild 






/ 




/ 






Moderate/severe 


0/40 


12/24 


N/A 


None/mild 






/ 




/ 






Not severe 


40/40 


18/24 


N/A 


Not severe 




/ 






x 




/ 


Severe 


0/40 


6/24 


N/A 


Not severe 




/ 






X 




/ 



CLJBW, creatinine clearance calculated using ideal body weight; CLM, renal clearance of D30G; CLP M15 , metabolic clearance of dapagliflozin to D30G; 
CLP olher , metabolic clearance of dapagliflozin to unmeasured metabolites; CLP ronal , renal clearance of unchanged dapagliflozin to urine; N/A, not applicable; QP, 
inter-compartmental clearance of dapagliflozin; V2M, central volume of distribution of D30G; V2P, central volume of distribution of dapagliflozin; V3P, peripheral 
volume of distribution of dapagliflozin; NA, not used during the model development. 

"Allometric scaling (BBWT,/70 kg) 3 ' 4 . "Creatinine clearance calculated using ideal body weight (IBW) was selected based on goodness-of-fit during develop- 
ment of the base model. "Normal, mild (Child-Pugh A), moderate (Child-Pugh B), severe (Child-Pugh C). "Including age, gender, smoking-status matched 
healthy controls. "Continuous covariates (i.e., age and CL crlBW ) were modeled as a proportional change model centered on the mean covariate (see supple- 
mental NONMEM code). 'For categorical covariates (i.e., sex and hepatic impairment) the indicator variable for the most common category was set to 1 and 
the proportional change for the less common category was estimated (see supplemental NONMEM code). 




15 30 60 90 



Baseline creatinine clearance (ideal body weight) (ml/min) 

Figure 2 Renal contribution to the metabolism of dapagliflozin 
to D30G in healthy subjects, and T2DM subjects with normal or 
impaired renal function. Solid squares, MB1 02007; open squares, 
MB1 02027; Line, population predictions. 

Model evaluation 

In the final model, all the structural parameters were esti- 
mated with good precision (% relative standard error < 20%). 
The addition of gender effect on CLP and CLM improved the 
precision in the estimate of the effects of baseline CL cr and 
severe HI (Child-Pugh Class C) on CLP M15 . 

The visual predictive checks confirmed that the model 
had good simulation properties (Supplementary Figure S3 



online). As seen in Figures 3 and 4, there was good agree- 
ment in model predicted and observed plasma and urine 
dapagliflozin and D30G concentrations for different catego- 
ries of renal function: CL < 30, 30-45, 45-60, 60-90, and 

cr ' ' ' 

>90 ml/min. However, given the high variability in urine con- 
centrations and the paucity of urine data in subjects with CL cr 
> 90 ml/min, the model predictions in this stratum should be 
interpreted with caution. In an external evaluation, the final 
model accurately predicted the plasma dapagliflozin and 
D30G concentrations in T2DM subjects with moderate renal 
impairment (Supplementary Figure S2 online). 

Model application 

Simulations with the final model were performed to gauge the 
effect of decreasing renal function on the exposure of dapa- 
gliflozin and D30G (Table 3). Compared with T2DM sub- 
jects with normal renal function (CL cr of 80-1 00 ml/min), the 
exposure of dapagliflozin (realized as AUC at steady state 
(AUC SS )) was 1 .26 and 1 .35 times higher than subjects with 
CL cr of 50-79 and 30-49 ml/min, respectively. The increase 
in D30G exposure exceeds the increase in dapagliflozin 
exposure (1 .32 times at CL cr of 50-79 ml/min and 1 .57 times 
at CL cr of 30-49 ml/min). Similar changes were seen when 
patient characteristics of the two phase I studies were used 
for simulations. 

DISCUSSION 
Main findings 

The semi-mechanistic model indicated that the influence of 
both renal and hepatic function on the CLP M15 was the most 
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Figure 3 Visual predictive check: plasma dapagliflozin and D30G by renal function. CL cr , creatinine clearance. Points, observed data; solid 
line, median of the observations; dotted lines, 5th and 95th percentiles of the observations; shaded areas, 95% confidence interval for 
simulated data (1 ,000 simulated data sets) for the corresponding percentiles. 




Figure 4 Visual predictive check: urine dapagliflozin and D30G by renal functions. Points, observed data; solid line, median of the 
observations; dotted lines, 5th and 95th percentiles of the observations; shaded areas, 95% confidence interval for simulated data (1 ,000 
simulated data sets). 



likely explanation of the altered dapagliflozin and D30G PK 
observed in the subjects with renal or HI. The fraction of 
D30G formed via the renal route decreased from 40-55% in 
subjects with normal renal function (CL cr >80ml/min) to -10% 
in subjects with severe renal insufficiency (CL cr = 13ml/min) 
(Figure 2).With severe HI, reduction of 42% in the CLP M15 
and 29% in the CLM resulted in D30G exposure comparable 
to that in healthy subjects. The subjects with moderate HI had 
nearly double the exposure to D30G. Therefore, differences 
in systemic exposure to D30G between the HI and healthy 
subject groups did not correlate directly with the severity of 
HI. The possible explanation for these findings is that D30G 
is mainly cleared via renal excretion, and subjects with severe 
HI may concomitantly present renal dysfunction. 12 



Effect of renal impairment 

Metabolism (primarily glucuronidation via UGT1A9) is the 
major elimination pathway of dapagliflozin: -60%, and <2% 
of an administered dapagliflozin dose is recovered from urine 
as D30G and dapagliflozin, respectively. 4 The tissue mRNA 
expression of UGT1 A9 in humans has been reported in the lit- 
erature. On the basis of the reverse transcriptase polymerase 
chain reaction data, the expression of UGT1A9 was found to 
be about eightfold higher in the human kidney as compared 
with that of the liver. 13 These data suggest the involvement of 
kidney in the metabolism of dapagliflozin. Our analysis sug- 
gested that the impact of renal impairment on the metabolism 
of dapagliflozin was manifested by the reduced metabolic 
clearance of dapagliflozin. Although there is evidence in the 
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Table 3 Effects of renal function on steady-state dapagliflozin and D30G exposure 



Renal 


Dapagliflozin AUC ss (ng*h/ml)° 






D30G AUC s! 


, (ng*h/ml) c 




function 


Phase l e 




Phase III' 




Phase I* 




Phase III' 




CL r " 


Median (95%CI) d 


Ratio" 


Median (95%CI) d 


Ratio" 


Median (95%CI) d 


Ratio" 


Median (95%CI) d 


Ratio" 


30-49 


71 1 (660, 777) 


1.26 


793 (644, 987) 


1.35 


1696 (1600, 1864) 


1.57 


2053 (1573, 2645) 


1.57 


50-79 


567 (485,647) 


1.07 


697 (563, 839) 


1.26 


1133(1012, 1317) 


1.36 


1283 (1010,1535) 


1.32 


80-1 00 


526 (260,822) 




515 (396,614) 




727 (319,1292) 




873 (664, 1118) 





AUC, area under the plasma concentration-time curve; CI, confidence interval. 

"Baseline creatinine clearance corrected for ideal body weight (ml/min). "Ratio of the median AUC compared with the median AUC for baseline creatinine 
clearance of 80-100ml/min. c Steady-state AUC. "95% CI for 100 simulated trials. °Covariates (baseline body weight, age, creatinine clearance corrected for 
ideal body weight and sex) from studies MB1 02007 and MB1 02027 (were used for simulations with a dose of dapagliflozin 1 0 mg daily. 'Covariates (baseline 
body weight, age, creatinine clearance corrected for ideal body weight and sex) from study MB1 02029 were used for simulations with a dose of dapagliflozin 
10mg daily. 



literature to suggest that chronic kidney disease contributes 
to decreased nonrenal metabolism, 14 the impact of chronic 
kidney disease on the liver metabolism of dapagliflozin is not 
known. Nevertheless, our simulation results along with the 
literature findings suggest that the kidney plays an important 
role in the metabolism of dapagliflozin via the UGT1 A9 path- 
way, and this may also explain the greater systemic exposure 
of dapagliflozin in subjects with moderate and severe renal 
impairment at steady state as compared with subjects with 
normal renal function. 

Effect of HI 

The model assigned the effects of HI on the formation 
(CLP M15 ) and clearance of D30G (CLM) during the stepwise 
covariate modeling. For technical reasons (linearized SCM), 
hepatic dysfunction was tested as two separate binary covari- 
ates (i.e., mild vs. moderate/severe, and mild/moderate vs. 
severe). The selection of mild/moderate vs. severe on CLP M15 
and mild vs. moderate/severe on CLM should be seen in con- 
text of the flexibility of the model, limited data, and the loss of 
information when categorizing HI. 15 Despite that all subjects 
contributed data on renal function, only six (9%) informed the 
model regarding moderate or severe HI. The results suggest 
that metabolic clearance of dapagliflozin to D30G (CLP M15 ) is 
preserved in patients with moderate HI. 

Absorption and oral bioavailability 

Misspecification of dapagliflozin absorption also resulted in 
misfit of D30G data during early sampling times. Therefore, 
considerable effort was made to characterize the absorption 
of dapagliflozin. A transit compartment model estimating both 
the mean transit time and number of transit compartments 
resulted in a flexible absorption model that well described 
the rich data during the absorption phase. The drawback was 
considerably longer computation time, especially when incor- 
porating the multiple-dose data. 

Absolute oral bioavailability and nonrenal clearance can 
be determined from modeling oral data when plasma and 
urine data are available in subjects with varying renal func- 
tion. 16-18 Estimating the absolute oral bioavailability, and 
thereby the absolute value of clearance, reduced the risk 
of biased estimates of apparent nonrenal clearance when 
only oral data were available. 19 The model-estimated abso- 
lute bioavailability of dapagliflozin was 82% (base model), 
84% (covariate model), and 86% (final model), which are 
in good agreement with the absolute oral bioavailability 



(geometric mean (90% CI): 78% (73-83%)) determined in 
healthy subjects following the concomitant administration of 
single intravenous 14 C-labeled microdose at time of maxi- 
mum concentration (T max ) of the extravascular (unlabeled) 
therapeutic dose. 3 

Model application 

Mild renal impairment is common in T2DM and presents a 
significant proportion of the population treated with dapa- 
gliflozin. Given the interplay in body weight, sex, and renal 
function, simulations were performed (using the covariates 
reported for each subject) to assess the effects of mild and 
moderate renal impairment on dapagliflozin and D30G 
exposure. As seen in Table 3, the simulations based on the 
subjects who participated in the phase I PK studies do not 
differ substantially from those for the subjects in the phase 
III study. The increases in D30G AUC SS of 1 .36 and 1 .32 for 
mild (50 < Cl cr < 80 ml/min), and 1.57 and 1.57 for moder- 
ate (30 < Cl cr < 50 ml/min) renal impairment, are compa- 
rable. However, the simulations using the characteristics of 
phase I subjects could have underpredicted the increase in 
dapagliflozin exposure compared with using the characteris- 
tics from the phase III study. Differences in body weights of 
subjects in the phase I and III studies should be considered 
when a model developed with mostly normal weight subjects 
is used to simulate into a new and obese population. 

In conclusion, the semi-mechanistic model developed to 
estimate the renal and nonrenal clearance of dapagliflozin 
presented a simplification of the interplay among biological 
systems; however, it provided a useful approach to evaluate 
the impact of renal and HI on the population PK of dapa- 
gliflozin. The results emphasized the role of renal metabolism 
of dapagliflozin to D30G, which is the major elimination path- 
way of dapagliflozin. The model estimated that the kidneys 
contributed 40-55% of the metabolism of dapagliflozin with 
normal renal function. In renal impairment and severe HI, 
both formation and elimination of D30G were decreased. In 
T2DM subjects, D30G systemic exposure (AUC ss ) increased 
more than dapagliflozin exposure. 

METHODS 
Study design 

Dapagliflozin and D30G plasma and urine concentra- 
tion-time data from a single- and multiple-dose (7 days of 
q.d. dosing) renal impairment study (NCT00554450), 78 
a 5-day single-dose HI study, 9 and a 52-week, phase III 
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(NCT00663260) clinical trial 20 in patients with T2DM and 
moderate renal impairment were available for analysis. 

The renal impairment PK/PD study evaluated the PK, PD, 
and safety and tolerability of oral dapagliflozin in eight healthy 
adult subjects (50 mg single dose) and in 32 T2DM subjects 
with normal, mild, moderate, and severe renal impairment 
(50 mg single dose followed by 20 mg daily dose for 7 days). 78 
The HI study compared the safety and serial PK of a single 
dose of dapagliflozin (10mg) in 18 subjects with HI and 6 
healthy subjects. 9 The phase 3 clinical trial investigated the 
glycemic efficacy, renal safety, and PK/PD of dapagliflozin in 
252 subjects with T2DM and moderate renal impairment. 20 

Each clinical trial was conducted in accordance to local 
and international ethics standards and was monitored for 
compliance with current good clinical practice. The appro- 
priate ethics and regulatory approvals were obtained before 
study initiation, and written informed consent was obtained 
from all participants. 

Analysis of dapagliflozin and D30G concentrations 

Analysis of dapagliflozin and D30G plasma and urine con- 
centrations was performed using high-pressure liquid chroma- 
tography with tandem mass spectrometry detection within the 
known period of stability. All determinations of dapagliflozin and 
D30G concentrations in human plasma and urine were gener- 
ated in analytical runs using appropriate calibration curves and 
quality control samples that met preestablished acceptance 
criteria. The standard curve ranges for dapagliflozin in plasma 
and urine were 1-500 and 1-1000ng/ml, respectively. The 
standard curve ranges for D30G in plasma and urine were 
1-500 and 10-5000 ng/ml, respectively. The between-run and 
within-run coefficients of variation for the analytical quality con- 
trols of both analytes in either matrix were <10%. 

Pharmacometric analysis 

Modeling approach The single/first-dose data were used 
to develop a semi-mechanistic (base) model to estimate the 
renal and nonrenal clearances of dapagliflozin. SCM building 
was used to evaluate the effects of renal function markers, HI, 
and other baseline covariates on the structural parameters 
(covariate model). After confirming that the single-dose models 
adequately described the multiple-dose data in subjects with 
T2DM and moderate renal impairment, the final model incorpo- 
rated important covariates identified in the previous analyses. 21 
The NONMEM code (ICON Development Solutions, Ellicott 
city, MD), example data set, and summary of data used in the 
model development and evaluation are provided in the Supple- 
mentary Material. 

The compartmental models were parameterized as clear- 
ance and volume terms using a parameterization that is sta- 
ble for different types of estimation methods. The structural 
model parameters were described as follows (Eqs. 1 and 2): 



MU,. = logP TN 



(MU, +il P ) 



(D 



(2) 



P^j is the typical value of parameter P, MU is a function of 



is a realization of a normally distributed ran- 
dom variable with zero mean and variance co 2 ,. Interindividual 
variability was reported as approximate percent coefficient of 
variation, and significant correlations between random effects 
(identified in the base model using nonparametric estimation) 
were retained in the subsequent models. 

Plasma and urine concentrations of dapagliflozin and 
D30G were analyzed simultaneously, but four separate 
residual error models were developed to describe the differ- 
ences between the individual predictions and the observed 
dapagliflozin and D30G concentrations in plasma or urine. 
The correlation in residual errors for dapagliflozin and D30G, 
when measured in either plasma or urine at the same time, 
was estimated (i.e., replicate residual error). 22 As an exam- 
ple, the residual error model for urine dapagliflozin is shown 
(Eqs. 3 and 4). 



W 



dapa^ 



IY 



'V-dapa^ ' ^prop.dapa,; 



^add^apa^ , 



(3) 



-■Y w 



dapa 0 



>a„ ' -\J^ ® replicate,, 



replicate^ c replicate L , 



(4) 



Y jjdapa , the /th dapagliflozin observation in urine in the /th 
individual^ was described by the corresponding model pre- 
diction ( V/,y,dapa u ), the standard deviation of dapagliflozin 
in urine (W dapau ), 6 proPidap% , 6 addidapV and the proportion of 
residual error common to both dapagliflozin and D30G in 
urine ( e replicat% ) • e ~ /V(0, 1) and £ ^(0,1) are 
residual random variables that are normally distributed with a 
mean of 0 and variance of 1 (fixed). is parameterized 

in terms of proportional (9 propdapau ) and additive (0 adddapau ) 
components. 

The first-order condition estimation method with tt-e inter- 
action in NONMEM (version 7.1.2) was used to obtain the 
parameter estimates. Standard errors were generated using 
Monte Carlo importance sampling EM assisted by mode a 
posteriori estimation. Xpose (version 4.3.2) (GNU Lesser 
Public License) was used to process the NONMEM output. 

Base model 

The absorption of dapagliflozin was described by a transit 
compartment model. 23 /f TRANSn - was calculated from the esti- 
mate of the mean transit time (i.e., the time for a drug mol- 
ecule to transit from the first to the last absorption site) and 
number of transit compartments (N + 1 ) (Eq. 5). The bioavail- 
ability was described using a logit model to ensure that the 
posterior individual estimates of bioavailability did not exceed 
100% (Eq.6). 

(5) 



MTT 



log 



TVBIO 



j Mrao,/ 



BIO 



' , ( TVBIO "I 

1+e io9 l™)b°-' 



(6) 



The availability of plasma and urine dapagliflozin and 
D30G supported three identifiable elimination pathways of 
dapagliflozin: (1) CLP renal , (2) CLP M15 , and (3) CLP olher (Eq. 7). 



P TV , P. is the value of parameter P for the /th individual, and CLP and CLM were p roportiona | t0 the C L cr and reported 



www.nature.com/psp 




Population Pharmacokinetics of Dapagliflozin in Renal or Hepatic Impairment 

van der Waft eta/ 



8 



as [L/77]/[mL/min[CL cr ]]. Metabolic clearances (i.e., CLP M15 and 
CLP other ) were scaled by the individual baseline body weight 
(WT.) using an allometric model, (wt,/70 kg), to account for 
the size of metabolizing organ(s). The FE and FM M15 were 
obtained from Eqs. 8 and 9. 



CLP, =CLP ,,+CLP M1c ,+CLP„, 



(7) 



The contribution of the kidney to dapagliflozin metabolism 
to D30G (FM M15rena| ) was the ratio of the renal component of 



ACLP.., 5 to the relative clearance ACLP.., ■ 

M15 M15 



FM, 



ACLP,. 



ACLP., 



(13) 



_ CLP renaU 



FE,, 



CLP, 



d-P M15/ 
CLP, 



(8) 



O) 



Covariate model building 

Covariate relationships with structural parameters were eval- 
uated visually and statistically using empirical Bayes esti- 
mates vs. covariates and stepwise general additive modeling. 
Only the covariate relationships significant in the stepwise 
general additive modeling were included in the SCM. The 
covariate-PK relationships that were evaluated are listed in 
Table 1 . The linearization SCM tool was used for forward and 
backward selection. 24 In this method, the individual predic- 
tions were linearized with respect to the e values (random 
effects). To implement the linearized SCM method, categori- 
cal covariates (e.g., Child-Pugh classification) were dichoto- 
mized (Table 2). 

The renal contribution to the metabolism of dapagliflozin to 
D30G (FM M15renal ) is calculated using Eqs. 10-13. 

The relative clearance for dapagliflozin to D30G (ACLP M15 ) 
for changes in baseline CL was calculated using the rela- 



tionship between CL and CLP, 

~ cr IV 

ate model: 



as estimated in the covari- 



Model selection and evaluation 

Improvements in goodness-of-fit were assessed by compar- 
ing the changes in the NONMEM objective function values 
(OFVs), improvements in diagnostic plots, reductions in 
interindividual variability of structural model parameters and 
residual error, and acceptable predictive performance. 

The simulation properties of the models were evaluated 
by visual predictive checks, a graphical assessment of the 
agreement between simulations from the model and the 
observations. To confirm that the model could predict plasma 
and urine dapagliflozin and D30G concentrations over the 
range of renal functions in the data, CL cr was used as the 
independent variable for visual predictive checks. To account 
for the different dose and covariate effects of dapagliflozin, 
the prediction and variance corrected visual predictive 
checks were used. 25 

The effects of impaired renal function on the steady-state 
exposure of dapagliflozin and D30G were realized by calcu- 
lating the AUC ss using Eqs. 14 and 15. The actual baseline 
body weight, baseline CL cr , and sex of the subjects in the 
phase I PK and III PK/PD studies were used in the simula- 
tions. In the simulations, a daily dose of dapagliflozin 10mg 
was administered to all subjects. 



Dapagliflozin AUC = 



Dose BIO, 
CLP 



(14) 



ACLP M15 = 1 + 9 CLc CLP M15 ■ (CL cr , - CL crmean ) 



(10) 



FM.„ t . Dose BIO 
Dapagliflozin 3 - O - glucuronide AUC =- M15J 



CLM, 



(15) 



where 0 CL CLP M15 is the proportional change in CLP M15 , 
CL . is the individual baseline CL„. and CL„ r _„ is the 

cr,; cr' cr,mean 

mean baseline CL cr . The nonrenal contribution to ACLP M15 
(ACLP ) was calculated using the estimate of 



0C|CLP M 15 



, was calculated using 

M15, nonrenal' ° 

from the final model (see Table 1) and the 
mean baseline CL of 80.14ml/min for the two clinical 

cr 

pharmacology studies based on the assumption that the 
kidneys do not contribute to the metabolism at zero level 
of CL „: 



ACLP, 



M15,nonrenal 



1 + e CLc CLP M15 (0-CL crmean ) = 0.598 (11) 
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The renal component of ACLP M15 (ACLP M15renal ) was the dif- 
ference between the ACLP M15 for the range of baseline CL cr/ 
values (Eq. 10) and the nonrenal effect (Eq. 11): 



ACLP. 



M 15, renal 



= ACLP U1S -ACLP U 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

y Dapagliflozin is a first-in-class, orally active 
inhibitor of the human renal sodium-glucose 
co-transporter 2 and is in development for the 
treatment of T2DM. Until now, no article has 
been published reporting the population PK of 
dapagliflozin and its major metabolite D30G. 



WHAT QUESTION THIS STUDY ADDRESSED? 

S This study evaluated the effects of renal or HI 
on the PK of dapagliflozin and D30G and also 
identified that both renal and HI decreased 
CLP M15 and CLM. 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

S The kidneys contributed 40-55% of dapagliflozin 
metabolism in subjects with normal renal 
function. This article is the first analysis using 
a semi-mechanistic model to determine the 
renal and hepatic contributions on dapagliflozin 
metabolism and presents a useful approach to 
evaluate the impact of renal and hepatic func- 
tion on drug PK. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

S The study results may promote model-based 
analyses and the consideration of both the ex- 
cretory and metabolic capacity of the kidney 
when investigating the effects of renal impair- 
ment on pharmacokinetics. 
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